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We analyze time-resolved spontaneous emission from excitons confined in self-assembled InAs 
quantum dots placed at various distances to a semiconductor-air interface. The modification of the 
local density of optical states due to the proximity of the interface enables unambiguous determina- 
tion of the radiative and non-radiative decay rates of the excitons. From measurements at various 
emission energies we obtain the frequency dependence of the radiative decay rate, which is only 
revealed due to the separation of the radiative and non-radiative parts. It contains detailed infor- 
mation about the dependence of the exciton wavefunction on quantum dot size. The experimental 
results are compared to the quantum optics theory of a solid state emitter in an inhomogeneous 
environment. Using this model, we extract the frequency dependence of the overlap between the 
electron and hole wavefunctions. We furthermore discuss three models of quantum dot strain and 
compare the measured wavefunction overlap to these models. The observed frequency dependence 
of the wavefunction overlap can be understood qualitatively in terms of the different compressibility 
of electrons and holes originating from their different effective masses. 

PACS numbers: 78.67.Hc, 42.50.Ct, 78.47.-p 
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I. INTRODUCTION 

Semiconductor quantum dots (QDs) are nanoscale 
solid state structures that provide three-dimensional 
quantum confinement of otherwise mobile charge carri- 
ers. Self-assembled QDs of InAs embedded in GaAs pro- 
vide confinement for both electrons and holes in a direct 
band gap semiconductor. Hence, they are optically active 
with the benefits of a high quantum efficiency and com- 
patibility with existing semiconductor technology. These 
properties make the QDs highly promising light sources 
for novel optical devices including optical quantum in- 
formation devices^. This has led to an increasing inter- 
est in the quantum-optical properties of QDs, and major 
achievements include the demonstration of the Purcell ef- 
fect for QDs in solid-state cavities^ and strong coupling 
between a single QD and the optical mode of a cavity-^. 
Very recently also electrical tuning of such quantum pho- 
tonics devices was demonstrated^, which is a significant 
milestone towards practical all-solid state cavity quan- 
tum electrodynamics devices. 

Despite the recent progress, a thorough understand- 
ing of the dynamics of light-matter coupling for QDs in 
nanostructured photonic media is still lacking. Such an 
understanding is required for quantitative comparisons 
between experiment and theory. The problem is two- 
fold, i.e. a detailed understanding of both the optical 
part and the electronic part is required. The optical part 
is described by the local density of optical states (LDOS) 
expressing the distribution of modes that the QD can ra- 
diate to, while the electronic part is determined by the 
exciton wavefunction for the QD. Here we will investigate 
an optical system where the LDOS can be calculated ex- 
actly, and use that to extract detailed information about 
the QD. Our experimental results are compared to a the- 



oretical QD model, and the effect of QD size, material 
composition, and strain is investigated. Such quantita- 
tive comparisons of experimental data to simple theoret- 
ical QD models are much needed in order to assess the 
full potential of QDs in nanostructured media for, e.g., 
single-photon sources^, low-threshold lasers 8 , or sponta- 
neous emission controlsiilii. 

When interpreting spontaneous emission decay curves 
from QDs, it is often implicitly assumed that the QDs 
primarily decay through radiative recombination, while 
non-radiative processes are negligible. Unfortunately, 
this assumption is not generally valid, and omnipresent 
non-radiative processes must be considered. Only few 
experiments have addressed this issue. Robert et al. es- 
tablished an upper bound on the contribution from non- 
radiative processes of 25% by measuring the ratio of the 
bi-exciton to exciton emission intensity at saturation 11 . 
Quantitative measurements of the radiative and non- 
radiative decay rates of QDs were only carried out re- 
cently using a modified LDOS both for colloidal QD a 12 i 13 
and for self- assembled QDs 14 . Precise measurements of 
the radiative decay rates are essential since nanophotonic 
devices rely on the ability to manipulate the radiative 
processes, while non-radiative recombination leads to loss 
in the system. 

As first pointed out by Purcell^ the radiative decay 
rate of an emitter is modified inside a structured dielec- 
tric medium, which is due to the modification of the 
LDOS. In early experiments by Drexhage 1 ^, this effect 
was experimentally demonstrated by positioning emitters 
in the proximity of a reflecting surface. We have recently 
employed the modified LDOS near a semiconductor-air 
interface as a spectroscopic tool to extract radiative and 
non-radiative decay rates and from that infer the over- 
lap between the electron and hole wavefunctions 1 ^. This 
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technique relies on the fact that the radiative decay rate 
is proportional to the LDOS, while the non-radiative de- 
cay rate is unaffected. In the present paper we expand on 
our previous work in particular by comparing the mea- 
sured radiative decay rate to theory, which requires a 
detailed model of the QD electron and hole wavefunc- 
tions. We have measured the radiative decay rate at 
different emission energies, which reveals the dependence 
of the QD optical properties on its size. We review the 
Wigner-Weisskopf theory of spontaneous emission from 
QDs, predicting an exponential decay of the exciton pop- 
ulation and the LDOS is derived for the applied interface 
geometry using a Green's function technique. We fur- 
thermore show that the radiative decay rate of a QD in 
a homogeneous medium is proportional to the square of 
the overlap between electron and hole wavefunctions and 
calculate the frequency dependence of this overlap us- 
ing a simple two-band model of the QD. The QD model 
is discussed in details and compared to our experimen- 
tal data employing realistic parameters as input to the 
theory. The pronounced size dependence of the electron- 
hole wavefunction overlap is found to originate from the 
differences in effective mass and binding energy of the 
electron and hole. Furthermore, we investigate three dif- 
ferent strain models for the QD and compare their pre- 
dictions to experiment thereby providing valuable insight 
on the complex strain mechanisms of self-assembled QDs. 

This paper is organized as follows: In section |TT] 
we present the experimental method and in section IIIII 
the experimental results. In section IIVI we discuss the 
Wigner-Weisskopf model for spontaneous emission and 
derive the relation between the radiative decay rate, the 
LDOS, and the wavefunction overlap. In section [V] we 
combine the analytical expressions for the radiative rate 
with the numerical results for the wavefunction overlap 
and compare theory with experiment. Finally, we present 
conclusions in section [VTl 



II. EXPERIMENTAL TECHNIQUE FOR 
DETERMINING THE RADIATIVE DECAY 
RATE OF QUANTUM DOTS 



Spontaneous emission of a photon from a QD occurs 
when an electron-hole pair (an exciton) recombines, as 
illustrated schematically in Fig. [TJa). As will be shown 
rigorously in section IIVI the QD radiative decay rate 
r r ad(r, oj, e p ) in a structured environment is proportional 
to the projected local density of optical states (LDOS) 
p(r, cj,e p ), where the projection is along the direction 
e p of the transition momentum matrix element, which 
corresponds to the orientation of the transition dipole 
moment. The LDOS is modified in an inhomogeneous di- 
electric medium due to optical reflections at interfaces. In 
emission experiments, the total decay rate is measured, 
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FIG. 1: (a) Schematic band diagram illustrating the sponta- 
neous emission process in a QD. An electron is excited op- 
tically from a valence band to a conduction band wetting 
layer (WL) state and the generated electron and hole relax 
to the lowest energy QD state on a picosecond timescale. 
The electron can subsequently decay by either radiative or 
non-radiative recombination with rates T ra d and T nra d re- 
spectively, (b) Schematic illustration of the sample under 
investigation. InAs QDs are embedded in GaAs and posi- 
tioned at different distances z to the GaAs-air interface, (c) 
The LDOS as a function of distance z to a GaAs-air interface 
for a dipole orientation parallel (solid curve) or perpendicular 
(dashed curve) to the interface. 



which can be expressed a; 



,17 



r(r,c,e p ) = r h r :f(u) p{r >") e ^ +r nrad ("), 



(1) 



where phom{u) is the density of optical states for a ho- 
mogeneous medium, and T nrat i(w) is the rate for non- 
radiative recombination, lu is the emission frequency and 
thus fkv the emission energy, and r the position of the 
QD. Non-radiative recombination is due to intrinsic QD 
processes and thus independent of the LDOS. TJ?°™(cj) is 
the radiative rate that the QD would exhibit in a homoge- 
nous medium without any boundaries. In our case, the 
refractive index of the medium is n = 3.5 corresponding 
to that of GaAs. Investigating T^°™(w) in detail pro- 
vides valuable insight into the properties of the exciton 
wavefunction confined in the QD potential. 

The exact nature of non-radiative recombination in 
QDs is not yet fully understood. It is often implicitly 
assumed that non-radiative recombination is negligible, 
but as we will see in the following even for very weak ex- 
citation intensities this is not a valid assumption. Possi- 
ble non-radiative processes include surface recombination 
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FIG. 2: (Color online) (a) Topographic atomic force micro- 
graph depicting a surface area of 1 x 1 (im of uncapped QDs 
on the unprocessed wafer. The color scale runs from run to 
20 nm. (b) Histogram of the QD height measured by analysis 
of the AFM data. The blue line is a fit to the histogram data 
using a log- normal distribution as discussed in the text. 



at the interfaces between the QD and the surrounding 
semiconductor material, Auger processes, and trapping 
of electron and/or holes at defects^, and any first princi- 
ples calculation of these effects is a tremendous task. Re- 
liable ways of extracting the radiative and non-radiative 
parts of the decay rate are therefore essential. 

The radiative and non-radiative decay rates can be 
separated by time-resolved spontaneous emission mea- 
surements if the QDs are placed in an environment with 
a known LDOS, cf. Eq. (fT]). A planar interface be- 
tween two regions with different refractive indices is the 
most simple example of such an inhomogeneous dielectric 
medium^. For this particular geometry, the LDOS can 
be calculated exactly and without any free parameters. 
Here we employ the interface between GaAs (n = 3.5) 
and air (n = 1) as illustrated in Fig. QJb). We calculate 
the LDOS by a Green's function technique and the re- 
sults for dipolc orientations parallel or perpendicular to 
the interface are shown in Fig. [ljc) . We stress that no 
assumptions need to be made about, e.g., the QD den- 
sity and the excitation beam profile in order to employ 
this experimental technique, as opposed to alternative 
ways of determining the radiative decay rate such as by 
absorption spectroscop} ^ 21 ' 22 . 



III. MEASUREMENTS OF SPONTANEOUS 
EMISSION DECAY RATES NEAR A 
SEMICONDUCTOR-AIR INTERFACE 

The starting point of our investigations is a GaAs wafer 
grown by molecular beam epitaxy (MBE) . The QDs were 
grown using the Stranski-Krastranow method on a (001) 
GaAs substrate. The growth sequence was 50 nm AlAs, 
610 nm GaAs, 2.0 monolayers (MLs) In As, 300 nm GaAs 
cap, and finally 2.0 ML In As. The AlAs layer was in- 
cluded for an optional epitaxial lift-off process, which 
was not employed here. Both In As layers formed self- 
assembled QDs, but only the embedded layer was op- 
tically active because non-radiative surface recombina- 
tion dominates for QDs at the surface. However, since 




FIG. 3: Illustration of the optical measurement setup. The 
sample is kept at 14 K in a cryostat and illuminated by a 
pulsed laser. The spontaneously emitted light is collected 
and can be directed either to a CCD camera for sample align- 
ment or to a spectrometer equipped with a fast single photon 
counting avalanche photodiode (APD) for time-resolved mea- 
surements. 
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FIG. 4: (Color online) (a) Measured fast decay rate at an 
emission energy of 1.204 eV versus the distance to the in- 
terface (black and red squares, upper trace). The blue line 
shows a theoretical fit to the data using the exact expression 
for the LDOS and Eq. fT} for a dipole orientation parallel to 
the semiconductor-air interface. The data points marked in 
red are excluded from the fit. The lower trace shows the slow 
decay rate which does not depend on the distance to the in- 
terface but retains a constant value given by the dashed red 
horizontal line. Inset: Decay curves measured for two identi- 
cal QD ensembles at two different distances to the interface 
at an emission energy of 1.204 eV. The data shown with 
blue squares (black triangles) are obtained for z = 109 nm 
(z = 170 nm). The red lines are fits to the data using a bi- 
exponential model and we extract the decay rates from these 
fits, (b) Values of the goodness-of-fit parameter Xr obtained 
for the individual bi-exponential fits. The values obtained are 
very close to the ideal value of unity. 



the QDs at the surface were fabricated under identical 
growth conditions as the embedded ones, the density of 
the active QDs can be determined by atomic force mi- 
croscopy (AFM) of the sample surface. Such an atomic 
force micrograph is shown in Fig. (Ha) . The image con- 
tains detailed information about the geometry and height 
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FIG. 5: (Color online) (a) Measured decay rate at an emis- 
sion energy of 1.204 eV versus the normalized LDOS (black 
and red squares) for a dipole orientation parallel to the 
semiconductor-air interface. The red line shows a linear fit 
where only the black points have been included. The red 
points have been omitted because at these distances from the 
interface a systematic deviation from theory was observed, as 
quantified by the data in the lower panel, (b) Correlation be- 
tween data and theory for the linear regression when system- 
atically excluding the points nearest the interface in the fit. 
The correlation parameter converges to a value close to unity 
(dashed blue line) when seven points have been excluded. 
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FIG. 6: (Color online) Decay rates as a function of distance 
to the interface for six different emission energies. Each curve 
has been vertically offset by 0.1 ns _1 for visual clarity. The fits 
have been obtained using a systematic exclusion of data points 
near the interface, cf. Fig. [5jb). The emission energies are 
1.170 eV (solid black squares), 1.187 eV (open black squares), 
1.204 eV (solid red circles), 1.216 eV (open red circles), 1.252 
eV (solid blue diamonds), and 1.272 eV (open blue diamonds). 
For all emission energies we note the excellent agreement with 
theory for z > 75 nm. 
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FIG. 7: (Color online) (a) The total decay rate (red squares), 
radiative decay rate (black triangles) and non-radiative decay 
rate (blue circles) as a function of QD emission energy. The 
radiative rate decreases with energy but the non-radiative rate 
simultaneously increases so that the total (measured) decay 
rate increases with increasing energy. The solid black line is 
the result of the theoretical model of the radiative decay rate 
using an aspect ratio of 1/2, omitting the wetting layer, and 
for an indium mole fraction of 0.95. The top scale shows the 
heights used in the calculation. The details of this calculation 
are presented in section IVl 



of the QDs at the surface, but the topography of the sur- 
face is convolved with the AFM tip shape function, and 
thus the width and exact geometry of the QDs cannot 
be extracted directly. The maximum height, however, is 
not subject to this effect and we used the AFM data in 
Fig. [5](a) to obtain a histogram of the heights recorded 
for 100 randomly selected QDs. The result is shown in 
Fig. HJb) . We fit the height histogram by a log- normal 

distribution given by f(h) = -^L^ cxp (- ( '" ( ^ A ' )2 ) 

with ho = 68, a = 0.32, and /i — 1.9, where ft, is a di- 
mensionless length scale normalized to 1 nm. We find an 
average aspect ratio defined as the height to lateral base 
diameter ratio of 1/3.6, but due to the convolution effect 
this is an upper bound. Since the transition energy of a 
QD depends sensitively on its height, the height distribu- 
tion function is central when testing theoretical models of 
the emission spectrum against experiments, as discussed 
in section [V] 

From the AFM data we found a QD density of 250 
u\m~ 2 , which corresponds to an average interdot distance 
of 60 nm. This number should be compared to typical 
length scales for various relevant QD interactions. Car- 
rier tunneling is negligible for distances beyond 15 nrn^ 
and the dipole-dipole interaction is only significant for 
distances close to the Forster radius, which is typically 
2-9 nmii. Therefore the measurements performed here 
provide ensemble averaged values of single QD properties 
with the advantage of an excellent signal-to-noise ratio in 
the measurements. 

The wafer was processed by standard UV lithography 
and wet chemical etching in five subsequent steps with 
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FIG. 8: (Color online) (a) The fast decay rate obtained from 
measurements on the unprocessed wafer (solid black squares) 
and the total decay rate in a homogenous medium (open red 
circles) obtained using the rigorous separation of radiative 
and non-radiative homogeneous decay rates. Evidently, the 
frequency dependence of the normalized LDOS (solid blue 
line) for the unprocessed wafer (z — 302 nm) results in faster 
decay rates for QDs on the unprocessed wafer than would be 
the case in a homogeneous medium, (b) The (normalized) 
inhomogeneously broadened emission spectrum of the QDs 
obtained under weak excitation. 



nominal etch depths of 160 nm, 80 nm, 40 nm, 20 nm, and 
10 nm by which we obtained 32 fields with specific dis- 
tances from the QDs to the semiconductor surface. The 
32 fields were nominally equidistantly spaced with 10 nm 
spacing. The wet etching was done using an etchant com- 
prised of H3PO4 (85%), H 2 2 (30%), and H 2 in the 
ratio 3 : 1 : 60, which has an etch rate on GaAs at room 
temperature of 1 nm/s. We found that this etchant re- 
sults in surfaces of good optical quality with low surface 
roughness. A schematic illustration of the resulting sam- 
ple is shown in Fig.QJb). Finally we measured the actual 
distance from the QDs to the semiconductor surface by 
using a combination of secondary ion mass spectroscopy 
and surface profiling from which we found typical depth 
uncertainties of ±3.0nm. 

The experimental setup is illustrated in Fig. 03 The 
sample was kept at 14K and irradiated by a mode-locked 
Thsapphirc laser emitting 300 fs pulses at 1.45 eV, which 
corresponds to excitation of the wetting layer states of 
the QD ensemble. The repetition rate was 80 MHz and 
we used an excitation intensity of 7 kW/cm 2 , which cor- 
responds to less than 0.1 excitons per QD generated in 
the wetting layer per pulse, i.e. only the QD ground 
state is populated. Excitation of the WL states is advan- 
tageous since the same density of excitons is generated 
independent of sample thickness, which would not be the 
case for excitation in the GaAs barrier since the samples 
have different thicknesses. The pump configuration is il- 
lustrated in Fig. [Ha). The spontaneous emission from 
the QD ensemble was collected and then dispersed by 



a monochromator with a spectral resolution of 2.6 meV 
from which it was directed onto a fast silicon avalanche 
photodiode (APD). 

Representative examples of two decay curves obtained 
for two different distances to the interface are shown in 
the inset of Fig. [DJa). We model our data with a bi- 
cxponential decay and typically find a goodness-of-fit pa- 
rameter Xr~ °f 1 • 2 , which is close to the ideal value of 
unity, thus confirming the validity of the model. The im- 
portant parameter extracted from the fit is the fast decay 
rate, which is due to recombination of bright excitons in 
the QD, i.e., it equals the total decay rate discussed in 
the previous section. 

We measured the decay curves for 32 nominally 
equidistantly spaced distances to the interface. We found 
that for the two distances closest to the interface there 
was no detectable spontaneous emission due to the very 
close proximity of the interface and/or damage by the 
etching process. In Fig. [Ja) we show the extracted fast 
decay rates (upper trace) versus the measured distance 
from the QDs to the interface for the remaining 30 sam- 
ples. As shown in Fig. |4jb) we obtain values of Xr near 
unity for the bi-exponential fits for all distances. The 
solid blue curve is a fit to the fast decay rates using the 
calculated LDOS for a dipole oriented parallel to the in- 
terface. When omitting the seven data points closest to 
the interface, as discussed in detail below, we find an ex- 
cellent agreement with theory. Since the slow rate (lower 
trace) does not depend on the distance to the interface 
and therefore does not depend on the LDOS it must be 
dominated by non-radiative decay. It is attributed as due 
to the recombination of dark excitons and the dynamics 
will be the subject of a future publication 2 ^,. In the re- 
mainder of this paper we will consider only the fast decay 
rate. 

When plotting the fast decay rates as a function of 
the LDOS normalized to the density of states of a ho- 
mogeneous medium, a linear dependence is expected, see 
Fig.[^a). However, close to the interface (z < 75nm) the 
measured decay rate is found to be systematically larger 
than expected by theory. This deviation could be due to 
enhanced recombination rates induced by, e.g., scattering 
or impurities at the semiconductor surface. In order to 
exclude these effects in our analysis of intrinsic QD prop- 
erties, we systematically excluded the data points closest 
to the interface in the fit. We found that the linear re- 
gression correlation parameter 26 obtained from the fit in 
Fig. 02a) saturated close to the ideal value of unity when 
excluding the seven innermost data points, see Fig. 02b). 
These points were consequently abandoned in the analy- 
sis. For distances z > 75nm we find an excellent agree- 
ment between theory and experiment, which allows reli- 
able extraction of the radiative and non-radiative rates 
of the QDs. 

The linear fit in Fig. 0H a ) i s based on Eq. (CQ) and 
contains two free parameters, namely the homogeneous 
radiative decay rate r^°™ and the non-radiative decay 
rate T nrad . We obtain T^T = 0.95 ± 0.03ns- 1 and 
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Fnrad = 0. 1 l±0.03ns~ 1 at an emission energy of 1.204eV. 
For reference, the theoretical curve for a dipole oriented 
perpendicular to the interface (i.e. parallel to the growth 
direction) is shown in Fig. [BJc). Clearly, this model 
cannot fit the experimental data, which confirms that 
the QDs are oriented in the plane perpendicular to the 
growth direction, which was previously established by ab- 
sorption measurements 2 ^. We measured the decay curves 
at six energies across the inhomogeneously broadened 
emission spectrum. When extracting the fast rate from 
the fits, we obtain the curves shown in Fig. [5] The dif- 
ferent emission frequencies shift the curves along the ab- 
scissa and more importantly the amplitude and ordinate 
offsets are changing, which corresponds to changes in r ra d 
and T nra d- This shows directly the frequency dependence 
of rj!°™ and T nra d, which will be discussed below. 

The total decay rate in a homogeneous medium was 
extracted from Fig. by the method outlined above and 
the result is shown in Fig. [2a). The total decay rate in- 
creases with increasing emission energy, which could sug- 
gest that the radiative rate increases with energy, as has 
been reported for colloidal QDs 28 . However, the opposite 
turns out to be true for self-assembled QDs. In this case, 
r^°™(cj) is found to decrease with increasing energy, and 
the overall increase in the total rate is due to the pro- 
nounced increase in T nrac i(ui) with emission energy. It 
should be stressed that such variations in the radiative 
rate can be assessed only because a modified LDOS is em- 
ployed allowing to separate radiative and non-radiative 
contributions. The striking energy dependence of the 
radiative rate can be explained as being due to the de- 
pendence of the electron and hole wavefunctions on the 
size of the QD, which will be discussed and analyzed in 
detail in section fVl 

Here we want to stress some potential pitfalls in the in- 
terpretation of the frequency dependence of spontaneous 
emission decay rates from QDs. As already pointed out, 
it is decisive to include the effects of non-radiative re- 
combination and implement a technique that allows to 
separate it from radiative recombination. Thus, an a 
priori assumption of negligible non-radiative recombina- 
tion would erroneously lead to the conclusion that the 
radiative rate increases with energy. Furthermore, when 
extracting quantitative data for the QD decay rates, one 
has to be aware of the possible influence of the presence 
of interfaces. In a homogeneous medium the LDOS is 
proportional to the emission energy squared, but this is 
not the case in proximity of interfaces. In Fig. [5£a) we 
compare the measured total decay rate versus emission 
energy for an unprocessed wafer and compare to the total 
rate that QDs in a homogeneous medium would exhibit. 
The latter has been obtained using the LDOS technique 
explained above and provides "undisturbed" QD prop- 
erties. For the unprocessed wafer, the distance to the 
interface is 302 nm, and for this distance the LDOS in- 
creases with increasing energy, cf. the solid blue line in 
Fig. [5Ja). This frequency dependence of the LDOS mod- 
ifies the measured emission rates, which should be taken 



into account, and is an example of the importance of 
considering nearby interfaces in quantitative assessment 
of QD properties. Experimentally this issue could be 
solved by growing a very thick capping layer on top of 
the QDs. 



IV. THEORETICAL DESCRIPTION OF 
SPONTANEOUS EMISSION FROM QUANTUM 
DOTS DUE TO INTERACTION WITH THE 
QUANTIZED ELECTROMAGNETIC FIELD 

In this section, we give a theoretical description of the 
radiative decay of excitons in QDs. For sufficiently small 
QDs, the energy difference between bound states in the 
QDs is much larger than the Coulomb energy and the 
effect of the Coulomb interaction on the internal exci- 
ton dynamics becomes negligible 29 . This means that the 
electron and hole comprising the exciton may be consid- 
ered independent, which is the strong confinement model. 
Furthermore, we employ here a two-band description of 
the QD including the effects of a wetting layer and strain, 
which is sufficient to capture the essential properties of 
QD o 30 i 31 . Our objective here is to explore the validity 
of this model by a thorough comparison to our measure- 
ments of the radiative decay rate and the emission spec- 
trum, which is carried out in section [V] Thorough explo- 
rations of even simple QD models are much needed since 
complete microscopically correct QD models are outside 
reach both due to the lack of experimental knowledge 
about exact atomic composition and computational com- 
plexity. 

Spontaneous emission occurs due to the interaction of 
the exciton with the continuum of vacuum modes. A rig- 
orous description of spontaneous emission requires a fully 
quantum description where both the radiation field and 
the exciton states in the QD are quantized. We employ 
here the Wigner-Weisskopf model of spontaneous emis- 
sion, which is valid when the LDOS varies slowly with 
frequency over the linewidth of the QD. This is an ex- 
cellent approximation for the dielectric structures inves- 
tigated here, but may break down for QDs in photonic 
crystal leading to intricate non-Markovian dynamics 32 . 
The radiative coupling strength is determined by the 
electron momentum matrix element and depends further- 
more sensitively on the overlap between the electron and 
hole envelope wavefunctions that in turn gives rise to the 
frequency dependence of the radiative decay rate. We de- 
rive here this frequency dependence, which will be com- 
pared in detail to the experimental data in section IVl 



A. Wigner-Weisskopf theory of spontaneous 
emission from quantum dots 

According to effective mass theory, the solution to the 
Schrodinger equation for an electron in a solid is given 
by ^„(r) = F„(r)M nj0 (r), where F n (r) is the envelope 
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FIG. 9: Schematic illustration of the QD geometry used to 
calculate the envelope wavefunctions. We consider a lens- 
shaped QD with lateral base diameter D and height h con- 
sisting of In c Gai_ c As on a wetting layer (WL) . The symmetry 
axis z is indicated along with the radial directions x,y. 



function, u n ,o(r) is the Bloch function evaluated at the 
band edge q = 0, and n G {c, v} denotes the conduction 
(c) or valence (v) band. The envelope function is the 
solution to the effective mass Schrodinger-like equation 
governed by the Hamiltonian^ 3 - 



H (t) =H kin>n + V n (r), 



(2) 



where i?km,n is the kinetic energy operator, and V n (r) is 
the band confinement potential. Here we consider lens- 
shaped QD geometries as shown in Fig. ©. For the 

conduction band the effective mass is isotropic, so that 
ft 2 i 

-^kin = -2^ V ' ;S7J7) v > where m is the elementary 
electron mass and m„(r) is the effective mass. For the 
valence band the anisotropy of the effective mass must be 
taken into account, which is discussed in Appendices lAl 
and [Bj 

For the valence band m v < and V n (r) < 0, which 
leads to negative eigenenergies. In the electron-hole 
representation 3 ^ we define the hole in the valence band as 
a particle with positive effective mass m/ 1 (r) = — m v {v) 
subject to a positive confinement potential V/j(r) = 
— V v (r) yielding positive eigenenergies. Clearly, the en- 
velope function remains the same in the new represen- 
tation, i.e. Fh(r) — F v (r). For III-V semiconductors 
the valence band is comprised of degenerate bands with 
different effective masses. However for QDs strain lifts 
this degeneracy and we may consider the valence band 
as a single band (the heavy hole band). We discuss de- 
tails of the band structure in the presence of strain in 
Appendix 

We describe the light-matter interaction by the Hamil- 
tonian H'(r, t) = — — p ■ A(r, t), where q is the elemen- 
tary charge, p is the momentum operator and A(r, t) 
is the vector potential of the quantized electromagnetic 
field. The latter is given by 35 - 36 

A(r, t) = £ (MrWe^ + A* (r) t e^*) . 

(3) 

Here /i = (k, s) is the combined wavevector k and po- 
larization index s £ {1,2}, is the optical angular fee- 
is a normalization constant with eo 



quency, e. 



denoting the vacuum permittivity, e M is a unit vector 
in the direction of the polarization s, A„(r) is the spa- 
tial field distribution function, and a M and are the 
field annihilation and creation operators, respectively. In 
a homogeneous medium the field distribution functions 

are given by plane waves (r) = ^= y , where e r = n 2 is 
the relative static permittivity of the material. We will 
be working in the Coulomb gauge in which the scalar 
potential of the electromagnetic field is zero and the di- 
vergence of the vector potential vanishes. 

We consider the initial state \i) — |c) ® |0) with 
c) = F e (r)\u c , ) corresponding to a QD with one elec- 
tron promoted to the conduction band and no photons 
in the radiation field. The final state relevant for sponta- 
neous emission is = \v) <g) ll^) with \v) = Fh(r)\u Vt o), 
where a photon is radiated to the mode fj, while the 
excited electron has decayed to the valence band. We 
have deliberately not written the spatial dependence of 
|c) and \v) for visual clarity. We can expand the in- 
teraction Hamiltonian by insertion of complete sets to 
obtain an expression containing the raising and lowering 
operators of the electronic system, which we define as 
<7 + = \c)(v\ and er_ = |v)(c|, respectively. It is conve- 
nient to change to the interaction picture, in which the 
time-evolution of the raising and lowering operators are 
given by cr+(t) = <j+e luJot and <7_(t) = er_e~ iWQt , where 
we have introduced the energy of the exciton transition 
hbjQ. Furthermore, assuming that the spatial distribution 
functions are slowly varying on the scale of the wavefunc- 
tions they can be evaluated at the center position ro of 
the QD, which is the dipole approximation. The interac- 
tion Hamiltonian in the interaction picture now reads 



H'(r ,t) 



- -^-(v\p\c)a. ■ -e^A* (ro)ate^^, 

TOO ^ Wp ^ * 

P * 

where A M = lu^ — lu and we have omitted the two terms 
proportional to e ±4A,it since they are rapidly oscillating 
as a function of time, which is the rotating wave approx- 
imation. 

The general state vector of the system can be expanded 

as 



(5) 



and inserting into the Schrodinger equation in the inter- 
action picture leads to the following equation of motion— 



dt 



CeW = -2^ IHp|c)l 
x f dt'c e {t') (°° d^ P^lM e-^-^-t') 



(0) 



where we have included an integration over a Dirac delta 
function in frequency, and assumed that the momentum 
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matrix element is constant within the linewidth of the 
QD. p(r ,uj,e p ) is the projected LDOS defined as 



p(r ,w,e p ) 



,| a |A p (r )| 2 (7) 



where e p is the unit vector specifying the direction of 
(u|p|c). This direction is determined by the Bloch ma- 
trix element as discussed below. Since p(ro, w, e p )/u> in 
Eq. ([6]) is slowly varying over the linewidth of the emitter 
so that it can be evaluated at the emission frequency lvq 
and taken outside the integral. In this case the QD pop- 



ulation decays exponentially |c e (t)| = 
with the radiative decay rate given by 



-r ra d(ro,^o,e p )t 



r rad (r 0) o;o,gp) = ^ ? 



V | Hp | c) |2P(ro."Q.fep). (8) 



This is the Wigner-Weisskopf result for spontaneous 
emission from solid-state emitters. It states that the ra- 
diative decay rate is proportional to the projected LDOS 
and the momentum matrix element. In the following sub- 
section, we discuss the evaluation of these two terms. In 
the experiment, the number of photons emitted per time 
is measured, which is given by 

-(r rad +r,„ ad )t 



N(t) = aT rad e 



(9) 



where additionally the rate for non-radiative recombina- 
tion has been added, and a is an overall scaling param- 
eter determined by the detection efficiency and the total 
number of photons recorded during the measurement pe- 
riod. 



B. Evaluation of the projected LDOS and the 
transition matrix element 



The projected LDOS can be calculated using a Green's 
function technique. In terms of the Green's tensor 
G(r, r',oj), the projected LDOS is given b y 17 ' 37 

2uj 



(ep ■Im(G(r,r,u)) -e F 



(10) 



where cq is the speed of light in vacuum. The LDOS 
is a classical electromagnetic quantity obtained by solv- 
ing Maxwell's equations. However, it enters the quan- 
tum optical theory of light-matter interaction where it 
describes the local density of vacuum modes that spon- 
taneous emission can occur to. For the particular case 
of a semiconductor-air interface as considered here, the 
Green's tensor is obtained by solving the following closed 
expressio n 17 i 38 



G(r, 



where 



M n = 



8nk 2 



2k 2 - k 2 p 





dk p -^(M Q 



M r ), (11) 



2k 2 - k 2 p 






2{k 2 -k 2 z ) 



(12) 



and 



(fcV - k z rP) 







(fcV - k 2 z rP) 




(2k 






2 rP) 



(13) 



Here k = |k|, where k = ik p ,k z ^k^) is the k-vector in 
cylindrical coordinates, z is the distance from the QD 
to the interface, and r s (r p ) is the Fresnel reflection co- 
efficient for s-polarized (p-polarized) lighti 7 -. The result 
of the calculation for a GaAs-air interface is shown in 
Fig. He) 

We now consider the transition matrix element. Using 
the fact that the momentum operator is a differential 
operator (p = —ihV), that the Bloch functions u c / v (r) 
are orthogonal, which we will describe below, and that 
the envelope functions are slowly varying on the scale of 
a lattice parameter—, we obtain 



(«|p|c)| 2 «|(F,|F e )| 2 |K |p|^ )| 2 
= \{F h \F e )f x ?Z?±^M, 



(14) 



This important result states that the transition matrix 
element is given by the product of the electron and hole 
wavefunction overlap and the squared Bloch matrix ele- 
ment |(it U! o|p| u c,o}| 2 - The magnitude of the Bloch matrix 
element is a material parameter that is expressed in terms 
of the Kane energy E p (c)r^. c is the indium mole fraction 
in the In c Gai_ c As alloy. The 

In the Kane mode l 18 ' 40 the valence band Bloch func- 
tions are written as linear combinations of the basis func- 
tions \u x ), \u v ), and \u z ) that carry the symmetry prop- 
erties of p-orbitals. The specific linear combination de- 
pends on the q-vector of the envelope function. QDs 
grown by the Stranski-Krastanov technique are typically 
flat structures placed on top of a wetting layer, and 
quantization along the growth direction (z) is therefore 
dominating 41 -. As a consequence we can set q x = q y = 0, 
which is exact in the limit of a quantum well, leading 
to q = |q|e z . In this case the heavy hole Bloch function 
can be written as \u V fi) = \uhh) — ^ (\ u x) ± A u y))- The 
conduction band Bloch function has s-symmetry, so the 
Bloch functions for the valence and conduction bands are 
orthogonal as was used above. Furthermore, this means 
that the matrix element (u„.o|p|it c ,o}| 2 is non-zero only 
for p x and p y from which we conclude that the dipolc 
axis of the QD is perpendicular to the growth direction 
in agreement with our experiment. 



C. Frequency dependence of the radiative decay 
rate of quantum dots 

We are now in a position to put together the results 
of the previous sections and calculate the frequency de- 
pendence of the spontaneous emission decay rates in a 
homogenous medium. By insertion of Eq. (JT]) in Eq. ([S]) 
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and using the plane wave expression for the field distri- 
bution functions we readily obtain 

itTM = 

ihp| C >i 2 - E i 6 p • ^\ 2 ^ - (15) 

The sum over all optical modes p, can be converted to 
an integration over all k-vectors where the dispersion re- 
lation for a homogeneous medium w M = k^co/n is used. 
The sum over polarizations yields a factor of 2. Using 
also Eq. (fTl)) we obtain the important relation 

r rad>) = R 3 ^(c^K^MlFeM)! 2 . (16) 

Here we have indicated explicitly that the envelope wave- 
functions depend on the emission energy since varying 
the QD size, and thereby the emission energy leads to 
modifications of the wavefunctions. This effect will be 
discussed in detail in section [V] Eq. (fl"6)) is the key re- 
sult used to interpret the experimental measurements of 
the radiative decay rate presented in this paper. It fur- 
thermore allows extracting an experimental value for the 
overlap between the electron and hole wavefunctions. 



Quantum dot height [nm] 

16 12 109 8 7 6 5 4 3 
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FIG. 10: (Color online) (a) Measured spontaneous emission 
spectrum (black curve, normalized) and calculated spectrum 
(red curve) using the height distribution of the QDs from 
Fig. HJb) and implementing strain model 1. The resulting 
parameters have been optimized to fit the spectrum, and we 
find an aspect ratio of 1/6 and an indium mole fraction in the 
QD of 39%. (b) Measured electron-hole wavefunction overlap 
(black squares) and the theoretical calculation (solid red line) 
using the same parameters as in (a). 



V. COMPARISON BETWEEN EXPERIMENT 
AND THEORY FOR THE ELECTRON-HOLE 
WAVEFUNCTION OVERLAP 

In this section we calculate the electron-hole wavefunc- 
tion overlap and the QD emission spectrum and compare 
to our experimental results. We investigate to what ex- 
tent the QD heights obtained from AFM measurements 
of uncapped QDs, cf. Fig. can be used as input to 
the models. Furthermore, we systematically test the 
model with experiment by varying parameters such as 
the indium mole fraction, the QD aspect ratio, the wet- 
ting layer thickness, and the applied strain model within 
physically realistic boundaries. The model nicely repro- 
duces the decrease of the electron- hole wavefunction over- 
lap with energy that we observed experimentally. How- 
ever, our investigations lead to the conclusion that fur- 
ther knowledge of QD composition and more involved QD 
models would be required in order to reach full quanti- 
tative agreement between experiment and theory. 

We describe the QD in the (p, z)-plane as one quarter 
of an ellipse with a fixed aspect ratio. Here p denotes 
the radial direction in cylindrical coordinates. The QD 
geometry is solved numerically in a large simulation area 
which spans 160 nm in the p-direction and 80 nm in the 
z-direction ensuring that the proximity of the boundaries 
have no effect on the results. Further details on the nu- 
merical procedure is provided in Appendix [B] We use 
the commercial finite element software package COM- 
SOL with an adaptive mesh to solve the effective mass 
equation (Eq. (|B3[) h For each QD height we obtain the 
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FIG. 11: (Color online) Calculated energy dependence of 
the squared electron and hole wavefunction overlap using the 
three different strain models discussed in the text. The QD 
aspect ratio is 1/2 and no wetting layer was included. The 
curves are calculated for strain model 1 with c = 0.51 (blue 
line), strain model 2 with c = 0.95 (red line) and strain model 
3 with c = 0.46 (black line). The data points show the ex- 
perimentally determined wavefunction overlap for these three 
indium mole fractions using the same color coding. Clearly, 
strain model 1 does not describe the experimental data well, 
while both models 2 and 3 lead to very good agreement. The 
inset shows the dependence of the QD emission energy on 
height for the three different models using the same color 
coding as in the main figure. 
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FIG. 12: (Color online) The effect of the thickness of the wet- 
ting layer on the calculated electron-hole wavefunction overlap 
when varying between ML and 6 ML in steps of 1 ML. Here 
we have omitted strain and otherwise used the same param- 
eters as in Fig. 1111 The experimental data for c = 0.46 are 
shown as black squares. 
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FIG. 13: (Color online) Color plots in the radial plane of the 
amplitude of the wavefunction for electrons and holes. The 
parameters corresponding to the calculation for strain model 
2 in Fig. [TT] were used, (a) Electron wavefunction for a QD 
height of 4 nm. (b) Hole wavefunction for a QD height of 
4 nm. (c) Electron wavefunction for a QD height of 8 nm. 
(d) Hole wavefunction for a QD height of 8 nm. 



transition energy and given the height distribution func- 
tion measured by AFM, c.f. Fig. [21 we calculate the 
emission spectrum. 

The fact that we observe transitions involving heavy 
holes motivates a further investigation of the analogy be- 
tween QDs and quantum wells. In particular the effect of 
strain on the electronic level structure is well-understood 
for quantum wells^. Strain has a significant effect on 
the level structure also for QDs, but is often discussed 
in qualitative terms only due to the mathematical com- 
plexity and lack of experimental input on the exact QD 



geometry and composition. Thus, we will in the follow- 
ing attempt to model the strain properties of the QD 
similar to the case of a quantum well and compare the 
results to experiment. In a quantum well of InGaAs in 
GaAs the compressive strain in the plane of the quantum 
well leads to an expansion in the direction perpendicu- 
lar to the plane. This biaxial strain can be decomposed 
into a hydrostatic and a shear component. However, as 
opposed to the case of a quantum well a QD cannot ex- 
pand freely in the growth direction which suggests that 
the hydrostatic component may dominate. Therefore, we 
compare three different strain models including: 1) hy- 
drostatic and shear strain, 2) only hydrostatic strain, and 
3) no strain. Strain modifies the band offsets and the va- 
lence band effective masses, which is discussed in further 
detail in Appendix [AJ 

Using the experimental data on both the QD height 
distribution function, the emission spectrum, and the 
frequency dependence of the wavefunction overlap, we 
explore experimentally realistic parameters in order to 
match our experimental data. The first approach is to 
use the measured QD height distribution of Fig. [2] and 
include both shear and hydrostatic strain. Using opti- 
mized parameters corresponding to an aspect ratio of 
1/6 and an indium mole fraction of 39% we find a very 
good agreement with the emission spectrum, as shown in 
Fig. [TUT a). The electron-hole wavefunction overlap ex- 
tracted from the radiative decay rate offers an additional 
test of this parameter set. In Fig. [TUf b) the experimen- 
tally determined wavefunction overlap is plotted along 
with the resulting theory 4 ^. The theory predicts correctly 
that the wavefunction overlap decreases (increases) with 
emission energy (QD size), and the mechanism behind 
this effect is discussed in detail below. However, a clear 
systematic deviation between theory and experiment is 
observed, which turns out to be the case for all three 
strain models provided that the model parameters are 
constrained to optimally reproduce the measured spec- 
trum. We conclude from this result that the measured 
height distribution of uncapped QDs is not very useful in 
determining the actual confinement volume of the over- 
grown QDs that the experiments are performed on. This 
is due to complex redistribution and intermixing pro- 
cesses of indium and gallium will occur during growth 
and subsequent regrowt h 43 ! 44 , which are likely to modify 
the QD confinement potential and strain significantly. 

Thus abandoning at this point the interpretation of 
the measured height distribution, which is closely related 
to the emission spectrum, we focus on the electron-hole 
wavefunction overlap. Keeping the aspect ratio fixed at 
a realistic value of 1/2 and excluding for the moment the 
wetting layer, we obtain the overlap shown in Fig. 1111 
Here the only free parameter is the mole fraction of in- 
dium c in the QD and the different strain models are com- 
pared. Note that the experimentally determined wave- 
function overlap depends on c, since it enters through 
the Kane energy, cf. Eq. (fPo)) . We find that includ- 
ing both shear and hydrostatic strain leads to a disagree- 
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ment between experiment and theory. This demonstrates 
that the strain model developed for quantum wells fails 
in the case of QDs, which is sometimes assumed in the 
literature^. Good agreement between experiment and 
theory can be obtained when including only hydrostatic 
strain or no strain at all for c = 0.95 and c = 0.46, 
respectively. Judging from experiments available in the 
literatur o 43 i 44 both these values of c are reasonable for 
overgrown QDs, so there is no support for favoring one 
of the two surviving strain models. We are led to the con- 
clusion that in particular shear strain is less significant for 
QDs compared to quantum wells, although further micro- 
scopic details of QD composition and geometry would be 
required for a further investigation of these issues. 

Figure rT2] investigates the effect of the wetting layer 
thickness on the wavefunction overlap. In this case strain 
is omitted, the aspect ratio is 1/2, and the indium mole 
fraction 0.46. We find that for wetting layer thicknesses 
below 4 monolayer (ML) the wavefunction overlap is only 
slightly modified in the emission energy range of interest 
to the experiment. For a very thick wetting layer (6 ML), 
the QD wavefunction is modified such that the mono- 
tonic decrease of the wavefunction overlap with energy 
observed for all other thicknesses does not apply. This 
behavior can be understood as follows: for very thick 
wetting layers and small QDs a significant part of both 
the electron and hole wavefunctions are expelled from the 
QD giving rise to quantum well-like wetting layer states 
that can have an increased mutual overlap. For the QD 
sample of the experiment the wetting layer was on the 
order of 2 ML. 

The above discussions point to a number of subtleties 
associated with a quantitative comparison between ex- 
periment and theory. These are mainly related to lack 
of knowledge about intrinsic properties of the QDs. No- 
tably, the monotonic decrease in the electron-hole wave- 
function overlap with emission energy is found to be a 
very general and robust result for a large range of dif- 
ferent parameters. The generality of this result can be 
understood in a simple physical picture and is related to 
the differences in the electron and heavy hole effective 
masses. Thus, for any indium mole fraction, we have 
from Eqs. (|A10|) . (|A11|) . and (|A12|) that m hh:b > m e ,b 
and mhh,x,Tnhh,z > w e . This in turn leads to a smaller 
Bohr oh = 4 ^ e o e '- h radius for holes than for electrons, 
i.e. the hole wavefunctions are more compressible than 
the electron wavefunctions. Increasing the emission en- 
ergy corresponds to decreasing the QD size. The large 
QDs emitting at small energies have a relatively large 
electron-hole wavefunction overlap. Decreasing QD size 
(i.e. increasing emission energy) compresses the electron 
and hole wavefunctions, and since the quantum confine- 
ment effect influences the electron wavefunction more 
than the hole wavefunction, the wavefunction overlap is 
decreased. This effect is illustrated in Fig. [13] where the 
calculated electron and hole wavefunctions for two dif- 
ferent QD sizes are plotted. It is clearly observed that a 
reduction in the QD height leads to a compression of the 



hole envelope wavefunction while the electron wavefunc- 
tion extends further into the surrounding GaAs barrier. 



VI. CONCLUSION 

We have presented time-resolved measurements of 
spontaneous emission from self-assembled QDs near a 
semiconductor-air interface. The interface leads to a 
modification of the LDOS, which can be calculated in 
an exact model without any adjustable parameters. The 
excellent agreement between theory and experiment en- 
ables separating radiative and non-radiative decay rates 
whereby they can be determined with unprecedented ac- 
curacy. We reviewed the theory behind the experiment 
by calculating the spontaneous emission radiative decay 
rate in a full quantum model in which spontaneous emis- 
sion is described by Wigner-Weisskopf theory. The ra- 
diative decay rate is proportional to the projected LDOS, 
which was derived using a Green's function formalism. 

From our measurements of the radiative decay rate at 
different emission energies we extracted the frequency de- 
pendence of the overlap of the electron and hole envelope 
wavefunctions. The experimental data were compared to 
theory by solving numerically the Schrodinger equation 
for a QD potential including the effects of shear and hy- 
drostatic strain. From this model the spontaneous emis- 
sion spectrum, which is inhomogeneously broadened due 
to the different sizes of QDs making up the ensemble, 
and the electron-hole wavefunction overlap were derived. 
An attempt to model the emission spectrum using the 
QD height distribution obtained by AFM on uncapped 
QDs was unsuccessful leading to the conclusion that this 
height distribution does not properly reflect the micro- 
scopic confinement potential of overgrown QDs. Regard- 
ing the frequency dependence of the electron-hole wave- 
function overlap, we found good agreement between ex- 
periment and theory with reasonable assumptions about 
QD size, geometry, strain, and wetting layer thickness 
assuming purely hydrostatic strain or no strain at all. 
In contrast, systematic deviations were found when in- 
cluding both shear and hydrostatic strain. Although the 
numerical model employed here is too simple to reflect 
the microscopic details of the QD geometry such as, e.g., 
indium-gallium intermixing, it reflects this simple phys- 
ical picture very well. A more detailed comparison be- 
tween experiment and theory is limited by the lack of 
experimental input on the QD structure at the atomic 
scale, which would be required to verify more sophisti- 
cated QD models. 

We finally discussed how the striking frequency depen- 
dence of the radiative decay rate, and consequently of the 
electron-hole wavefunction overlap, can be understood 
in terms of a very simple physical picture: The heavy 
hole wavefunction is more compressible than the elec- 
tron wavefunction due to the larger effective mass. A 
reduction in the QD size therefore leads to a further lo- 
calization of the heavy hole wavefunction, while the elec- 
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tron wavefunction delocalizes into the surrounding bar- 
rier. This leads to a decrease in wavefunction overlap for 
increasing emission energy. 

Quantitative comparison of the experimental data to 
theory was limited by the lack of detailed experimen- 
tal input about the microscopic composition of the QDs. 
Combining the detailed optical experiments presented 
here with techniques to extract local material proper- 
ties of QDs, e.g., by high resolution transmission elec- 
tron microscopy (TEM), will be a very exciting future 
research direction that also will pinpoint the need for 
more involved theoretical models of the QDs. We believe 
that the technique presented here to directly access the 
light-matter coupling strength will have important appli- 



cations regarding proper design and characterization of 
solid state quantum photonic devices. 
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FIG. 14: Schematic illustration of the effect of strain on the 
band energies at q = 0. (a) No strain and (b) biaxial strain. 
Biaxial strain shifts both the conduction band and the heavy 
hole valence band closer to the continuum. The net effect is 
an increase in the band gap and therefore in the transition 
energy. Importantly, the degeneracy of the light and heavy 
hole states is lifted and the lowest energy transition involves 
heavy holes only. 



37 H. T. Dung, L. Knoll, and D.-G. Welsch, Physical Review 
A 62, 053804 (2000). 

38 M. Paulus, P. Gay-Balmaz, and O. J. F. Martin, Phys. 
Rev. E 62, 5797 (2000). 

39 L. A. Coldren and S. W. Corzine, Diode lasers and pho- 
tonic integrated circuits (Wiley, 1995). 

40 J. H. Davies, The physics of low- dimensional semiconduc- 
tors: An introduction (Cambridge University Press, 1998). 

41 A. D. Andreev and E. P. O'Reilly, Applied Physics Letters 
87, 213106 (2005). 

42 Note that in Fig.3.B of Ref. [14] the theoretical 
curve plotted is mistakenly \{Fh{ui)\F e {uj))\ rather than 
|(F,H|F e H)| 2 . 

43 I. Kegel, T. H. Metzger, A. Lorke, J. Peisl, J. Stangl, 
G. Bauer, J. M. Garcia, and P. M. Petroff, Physical Review 
Letters 85, 1694 (2000). 

44 N. Liu, J. Tersoff, O. Baklenov, A. L. Holmes, Jr., and 
C. K. Shih, Physical Review Letters 84, 334 (2000). 

45 O. Stier, M. Grundmann, and D. Bimberg, Physical Re- 
view B 59, 5688 (1999). 

46 I. Vurgaftman, J. R. Meyer, and L. R. Ram-Mohan, Jour- 
nal of Applied Physics 89, 5815 (2001). 

47 W. W. Chow and S. W. Koch, Semiconductor-laser funda- 
mentals (Springer, 1999). 



APPENDIX A: INFLUENCE OF STRAIN ON 
EXCITONS CONFINED IN QUANTUM DOTS 

Strain due to the lattice constant mismatch between 
InAs and GaAs is responsible for the formation of QDs 
during MBE growth in the Stranski-Krastranov growth 
mode. This means that the QDs are highly strained and 
this has significant impact on the electronic band struc- 
ture. The interplay between geometry, chemical compo- 
sition, and strain is complicated for QDs. During the 
growth, diffusion of In and Ga takes place, so that the 
resulting QD will consist of a significant fraction of Ga 



even if it is grown by pure InAs 43 . Furthermore, it has 
been reported that In is mainly concentrated in an in- 
verted cone inside the QD giving rise to a complex strain 
profile 44 . Complete knowledge about such complex de- 
tails is still lacking, and the purpose here is to introduce 
simple strain models and judge their validity by com- 
paring to experimental data. We consider a lens-shaped 
QD with a lateral extension, which is larger than the 
extension in the growth direction. For this reason we ap- 
proximate the strain of a QD by the model used in the 
case of a quantum well. This is further motivated by 
the fact that the QD is placed on top of a wetting layer, 
c.f. Fig. [9] The strain modifies the band offsets and en- 
ergy gap of the QD, and the valence band degeneracy 
is lifted so that the transition with the lowest energy in- 
volves only heavy holes. This is illustrated in Fig. ITU We 
will assume that the bulk GaAs surrounding the QD is 
unstrained, however only heavy hole bands are included 
here as well since the influence of band mixing is minor 
in the barrier where the electron and hole wavefunctions 
are strongly damped. 

One effect of strain is to shift the conduction and va- 
lence bands in energy. The strain describes the compres- 
sion or expansion of the crystal lattice, which in general 
is described by a tensor e nm . For biaxial strain, which 
describes the strain at planar heterojunctions, only the 
diagonal elements are relevant. We consider a thin layer 
of In c Gai_ c As with lattice constant oqd embedded in a 
GaAs barrier with lattice constant a,b- We assume that all 
strain is incorporated in the In c Gai_ c As layer and since 
a QD > &b the strain will be compressive. We havei^ 



~yy 



2C 



12 



c 



11 



(Al) 



where C\\ and C12 are the elastic stiffness constants (ma- 
trix elements of the stiffness tensor). These strain com- 
ponents lead to a change in band structure and thus a 
modification of band energies and effective masses, as de- 
scribed by the Pikus-Bir strain models. Biaxial strain 
can be decomposed into components of hydrostatic and 
shear strain. The hydrostatic compressive strain of the 
QD leads to a decrease of the band offsets exactly as for 
any hydrostatic compressive strain resulting from, e.g., 
a decrease of the temperature. For a thin strained epi- 
taxial layer, the energy can be lowered by compensating 
the in-plane compressive strain by an expansion in the 
z-direction (shear strain) . We obtain the following heavy 
hole valence/conduction band offsets^ 



AE hh 
AE r 



AE' V - 



where 



Re — ^c{^xx €yy ~t~ ^zz)i 



(A2) 
(A3) 

(A4) 

(A5) 
(A6) 



14 



TABLE I: Material parameters for In c Gai_ 
and valence band parameters respectively. 



; As at cryogenic temperatures used in this work. CB and VB indicate conduction 



Quantity 




Value for In c Gai_ c As 


Unit 


Refe: 


Lattice constant 


a 


5. 


6503 + 0.4050c 


A 


^5 


Band gap 


E g 


1.515 


- 1.580c + 0.475c 2 


eV 


45 


CB effective mass 


m e 


0.0667 - 


- 0.0419c - 0.00254c 2 


mo 


45 


Luttinger parameter 


7i 


1/[(1 


-c)/6.98 + c/20.0] 




45.46 


Luttinger parameter 


72 


1/1(1 


- c)/2.06 + c/8.5] 




45.46 


Luttinger parameter 


73 


VI(1 


- c)/2.93 + c/9.2] 




45,46 


CB Hydrostatic def. pot. 


o c 




8.013 + 2.933c 


eV 


45 


VB Hydrostatic def. pot. 


a v 




1.824 + 0.024c 


eV 


45 


VB Shear def. pot. 


b v 




-2.0 + 0.2c 


eV 


46 


Elastic stiffness constant 


C11 




1221 - 388.1c 


GPa 


46 


Elastic stiffness constant 


C12 




566 - 113.4c 


GPa 


46 


Static dielectric constant 






13.18 + 1.42c 




45 


Kane Energy 


E p (c) 




28.8 - 7.3c 


eV 


39.46 



and we have introduced a number of quantities denned 
in Table Q] AE' C (AE' V ) is the unstrained conduction 
(valence) band offset which constitute 60% (40%) of the 
band gap difference between bulk and QD, so that, e.g., 
AE' C = 0.6(^(0) - E' g (c)). The band gap of the strained 
QD is given by 



E g (c) = E'Jc)+P e + Q e + R e 



(A7) 



where E' g {c) is the unstrained band gap of the mate- 
rial. These effects on the band structure are illustrated 
in Fig. M 

Another important consequence of strain is that the 
effective heavy hole mass becomes highly anisotropic. In 
contrast the effective electron mass is not modified con- 
siderably since the conduction band is much more ener- 
getically isolated than the valence band o 18 i 46 . In our ex- 
periments the growth direction (z) is parallel to the [001] 
crystal axis, which is perpendicular to the wetting layer 
plane (x,y) oriented along [110] and [110] directions. We 



will use 



f [ooi; 



and 



i x = = m [110] 



Fur- 



thermore, since the crystal structure along the [110] and 
[110] axes are identical apart from a series of rotations, 
we have that m x = m y . 

The heavy hole masses for unstrained bulk In c Gai_ c As 
are given by 



hh,\\,b 



L hhA.b 



2 7l - 72 - 373 ' 

1 



7i - 2 72 



(A8) 
(A9) 



where 71, 72 and 73 are Luttinger parameters^. These 
are listed along with all other relevant QD material pa- 
rameters for this work in Table HI For In c Gai_ c As, 
we have 72 ~ 73 allowing for the axial approxima- 
tion. Here 72 and 73 are replaced by their average value 
7 = (72 + 73)/2, leading to an isotropic heavy hole mass 
valid for unstrained bulk semiconductors^ 



l hh,\\,b — m hh,A_,b 



1 



7i - 2 7 



(A10) 



These bulk effective masses are used to describe the GaAs 
barrier material surrounding the QDs. The strained 
heavy hole effective masses in the directions parallel with 
and perpendicular to the wetting layer plane are given 



m hh,\\,QD 



TUhh, _L,QD 



1 



71+7 
1 

71 - 27 ' 



(AH) 
(A12) 



showing that the parallel component is strongly modified 
by strain. 



APPENDIX B: NUMERICAL MODELING OF 
ENVELOPE WAVEFUNCTIONS 

We solve the effective mass equation for electrons in the 
conduction band and holes in the valence band in order 
to calculate the overlap of the envelope wavefunctions. 
The effective mass equations describing the electron and 
the hole have the same form, but the effective hole mass 
anisotropy must be taken into account. We therefore 
consider the anisotropic valence band problem which has 
the isotropic conduction band problem as a special case. 
For both electrons and holes the effective mass depends 
on position due to the different effective masses in the 
QD and in the surrounding crystal matrix. 

Using cylindrical coordinates (p, 4>, z) and the axial ap- 
proximation, the kinetic term in Eq. ([2]) reads 



H 



kin 



1 d 



P 



d 



2m \pdp\ m„|| (p, z) dp 
d_( 1 d 

dz \m n ^(p, z) dz 



1 d 2 



w„||(p, z)p 2 def) 2 



(Bl) 



Using separation of variables F(r) — R(p, z)<&((f) the ef- 
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fective mass Schrodinger equation can be written as 



-1 8 2 



c) 



R{p, z) dp 

m n \\(p,z)p 2 d 



R{p,z) dp \m n \\{p,z) dp 
m n \\(p,z)p 2 d 



-R(p,z] 



d 



+ 



R(p,z) dz \m nJ _(p, z) dz 
2m m n \\(p,z)p 2 



R{p,z) 



K 2 



-(V(p,z)-E). 



(B2) 



The left and right hand sides of this equation are inde- 
pendent and they must therefore equal a constant, i.e. 

^7^(0) — — I 2 - The solution of this equation is 
$(</>) = cicos(lcj>) + C2sin(l(f>) which by the boundary con- 
dition $(0) = <&(27r) implies that I must be an integer. 
We are considering the ground state transition and there- 
fore take I = 0. This leaves an equation describing the 
electronic motion in the {p, z)-plane. 



Eq. (|B2[) is solved numerically after being reduced to 
a dimensionless form in order to avoid numerical issues 
related to the very small factors (oc h 2 ) appearing in 
this equation. We define the new dimensionless quan- 
tities through p = k p p, z — k z z, R(p,z) = kftR(p, z), 
V{p,z) = V V(p,z), E = V E, and V = h 2 /(2m k 2 ), 
and furthermore take k p = k z = 1 nm so that all spatial 
dimensions are measured in units of nanometer. By this 
transformation we obtain 



1 d „, d ( 1 d -,_ 
-R(p,z) + — — -——R(p lZi 



m\\(p, z)pdp 



V(p,, 



I 2 



dp \m||(/5, z) dp 

d ( 1 d - 

-t-A t^jtzR&z) ) (B3) 

dz \mj_{p, z) dz ' 



R(p,S) = ER(p,z), 



my (p, z)p 2 



which is solved numerically using a finite clement 
method. 



